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Abstract 

The neutrinos detected from the next Galactic core-collapse supernova will contain valuable 
information on the internal dynamics of the explosion. One mechanism leading to a temporal 
evolution of the neutrino signal is the variation of the induced neutrino flavor mixing driven by 
changes in the density profile. With one and two dimensional hydrodynamical simulations we 
identify the behavior and properties of prominent features of the explosion. Using these results we 
demonstrate the time variation of the neutrino crossing probabilities due to changes in the MSW 
neutrino transformations as the star explodes by using the S-matrix - Monte Carlo - approach to 
neutrino propagation. After adopting spectra for the neutrinos emitted from the proto-neutron star 
we calculate for a Galactic supernova the evolution of the positron spectra within a water Cerenkov 
detector and the ratio of charged current to neutral current event rates for a heavy water - SNO 
like - detector and find that these detector signals are feasible probes of a number of explosion 
features. 

PACS numbers: 14.60.Pq, 97.60.Bw 
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I. INTRODUCTION 



The importance of neutrinos in the explosions of massive stars has long been recognized. 
This significance is coupled with their ability to carry information to us about the processes 
and conditions in the core of the supernova so that the neutrino signal from the next Galactic 
supernova will provide us with an opportunity to test the core collapse paradigm. The 
explosion begins when the runaway process of electron capture in the core a massive star 
leads to a rapid compression and collapse that is only halted when the degenerate pressure 
of the neutrons and thermal pressure of non-degenerate particles kicks in at super-nuclear 
densities. The rapid neutronization during the collapse leads to a large burst of neutrinos 
that identifies the beginning of the event. As nuclear densities are reached the mean free 
path of the neutrinos becomes shorter than the size of the proto-neutron star and the 
neutrinos become trapped. Thermal processes within the core create a thermal bath of 
neutrino-antineutrino pairs of all flavors. The neutrinos slowly diffuse from the core over a 
period of order 10 seconds, carrying 99% of the gravitational binding energy of the core with 
them. With such a large neutrino luminosity even a small number of neutrino interactions 
above the core can create an important impact on the explosion physics. Detection of the 
neutrinos from the next nearby supernova will offer us the opportunity to examine the 
internal evolution of the explosion since the neutrinos allow us to see down all the way to 
the proto-neutron star. For example, if the proto-neutron star collapses to a black hole, the 
resulting neutrino signal will be altered Q] • 

In addition, although these neutrinos are only weakly interacting they can have significant 
effects upon the nucleosynthesis that occurs in the supernova: in the outer layers of the 
star neutrinos can cause a transformation of the elements synthesized during the preceding 
thermonuclear burning period in the "neutrino process" 3|, it is thought that a neutrino 
driven wind occurs at late time in the supernova which may create sufficiently neutron rich 
material to produce the r-process elements [sl, and long duration gamma ray bursts (which 
are thought to be a rare type of supernova) produce elements in conditions governed by 



strong neutrino fluxes [y, [Tj. The results o: 



re tvp 



sensitive to the neutrino spectra ,lio|,liil. 



all these nucleosynthesis processes are quite 



For reasons pertaining to both hydrodynamics and to element synthesis it is important 
to understand the details of neutrino propagation in the supernova and the mixing between 
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neutrino flavors. Neutrino oscillations also alter what we observe but, simultaneous^ 

ai: 



ow the possibly of elucidating yet unknown neutrino oscillation physics, e.g. 
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also 
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161 ]. Supernova neutrino flavor transformation is a rapidly developing field with the 
neutrino background terms the subject of intense study 13, Q, Q- Realistic possibilities 
include a complete mixing of all flavors 17] or a partial oscillation between flavors isl. An 
important part of any future observation of supernova neutrinos will be to understand this 
)hysics. In parallel, phase effects due to multiple resonances have recently been recognized 



2A m d.stortmg the s,mple p.cture of „eutn,jopropagatio„ through the outer layers of the 

supernova that has been used in the past [12, Il3|, The explosion of the star leads to an 
evolution of the flavor mixing so that to appreciate the effects upon neutrino propagation it 
is important to understand the density profile. 

Observations of high pulsar velocities and polarized supernova light suggest an inherent 
asymmetry in the explosion mechanism. Though many possible mechanisms for generating 
;he asymmetry have been proposed, recent work by Blondin, Mezzacappa & DeMarino 
22] has identified an instability of the standing accretion shock that leads to large dipole 
and quadrupole moments. More recent work by Blondin & Mezzacappa 23| and Ohnishi, 



Kotake, and Yamada 



241 ] as well as Scheck et al. [25|] has both confirmed and furthered 



the understanding of this instability. An alternative mechanism, acoustic heating, has also 



been discussed by Burrows et al. 



26| . The temporal evolution of the neutrino flavor mixing 



and the subsequent variation of the neutrino signal means that it may be possible to to 
detect prominent features of the supernova density profile and to use this information to 
learn important information about the explosion. This possibility was demonstrated first by 



Schirato and Fuller 



271 ] who used a time- dependent, one- dimensional 1/r density profile to 



demonstrate that the supernova's forward shock wave reaches and disrupts the 'H' resonance 
transformation layer a few seconds after the core bounce leading to a detectable change in 
the charged current neutrino signal. This was followed and elaborated upon by Dighe & 



Smirnov 



Takahashi et al. 



and Fogli et al. [29!]. More recently Tomas et al. 
used a more sophisticated numerical model of a progenitor star to investigate the effect of 
both forward and reverse shocks on the neutrino signal. Tomas et al. also present one 
second of a two-dimensional simulation result to show how the density profile can be greatly 
complicated by a deformed forward shock followed by strong convection currents. But due 
to the intensive computational burden of supernova models it is numerically very taxing to 
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to watch these the explosion propagate outward until they affect neutrino transformations. 

In this paper we present calculations which link hydrodynamical simulations of time de- 
pendent density profiles and phase-retaining neutrino oscillations. Our presentation begins 
in section ini with the results from simulations of supernovae in both one and two dimensions 
where we artificially heat a density profile constructed so as to mimic the state of the su- 
pernova at the point where the accretion shock has stalled. Our two-dimensional explosions 
are heated aspherically so as to create an aspherical supernova. Though this 2D model lacks 
some structural features present in more sophisticated models that generate asphericity via 
standing accretion shock instabilities, we still obtain results that give a good idea of how 
density isotropy affects neutrino flavor transformations. In section III we pass a spectrum of 
neutrinos through the simulation results and show how the crossing probability is affected 
by the evolution of the profiles. We finish in section IV by presenting our calculation of 
two neutrino detector signals - the positron spectra within a water Cerenkov detector, and 
the ratio of charged current to neutral current event rates for a heavy water detector - that 
demonstrate our ability to extract information about the supernova. 

II. THE SN PROFILE 
A. The SN Simulations 

Our first task in determining the neutrino signal is to simulate the supernova so as 
to obtain the time evolution of the density profile. A supernova simulation is a complex 
and computationally intensive problem. At the present time there is not yet a robust, self- 
consistent model for core-collapse supernovae, and the most sophisticated multi-dimensional 
simulations take a very long time to compute. But we do not need to simulate the actual 
core collapse and formation of the proto-neutron star, these events occur deep within the 
core at a radius of 10^ cm or so, whereas the first neutrino resonance occurs at around 
10^ cm for neutrino energies of order 10 MeV, and the second is even further out. What is 
occurring in the core is irrelevant for our purposes, all we need from the simulation of the 
supernova is the propagation of the shock through the star. For this reason we employ a 
simplistic and artificially energetic numerical simulation to create explosions with a range of 
features that might be present in reahty. We do not intend this model to be realistic itself, 
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but rather a tool to help with the difficult problem of reconstructing density features of an 
exploding star based on changing neutrino signals detected at Earth. 



We use VH-1 jsil, a hydrodynamic code based on Woodward and Colella's 32| piecewise 
parabolic method. For both the one-dimensional and two-dimensional simulations we employ 
an exponentially coarsening radial grid that ranges from r = 10 km to r = 250, 000 km where 
r is the radius from the center of the star, and for the two-dimensional simulations an angular 
grid that ranges from 6* = to = tt where 9 is the polar angle. The innermost radial grid 
increment is set to 5r = 202.6 m and increases by 1.3% thereafter. By r = 10'^ km the 
radial increment has grown to 6r ~ 6.5 km and by r = 10^ km we reach 6r ~ 650 km. 
We map into the code a spherically symmetric progenitor proffie intended to represent the 
state of the supernova about 100 ms after the collapse. Beyond r = 200 km t he p rofile 
is that of a low metallicity, 13.2 Mq progenitor model developed by Heger et al |33| that 
we have allowed to collapse further within VH-1 so that the infall velocities approached 
free fall. Interior to r = 200 km we splice in a slow, outward-moving standing accretion 



shock proffie provided by Heywood 3^ that is similar to that found in 2^. Finally, inside 
r = 20 km we insert a dense, 7 = 2.5, polytropic core containing approximately 3Mq. The 
inner boundary condition inside the dense core at 10 km is reflecting. Figure ([1]) shows the 
full initial density profile. To drive the explosion, we mimic neutrino heating by inserting 
energy into the region above a = 100 km gain radius. The energy deposition rate per 
unit volume, dQ/dtdV, is proportional to the density of material, falls as 1/r^ and decreases 
exponentially with time t over a timescale r set to r = 0.5 s. For the ID simulations then, 

dQ _ p fr-rg\ ./^ 



The additional factors of {r — rg)/r and t are inserted so as to ramp up the energy deposition 
over both distance and time thereby avoiding the violent disruption that occurs if we deposit 
energy either too suddenly or within a small volume. For the 2D models we introduce an 
angular dependence into the energy deposition prescription so as to match the observation 
by Blondin et al. that small perturbations in standing accretion shock models can 
ultimately lead to aspherical shock modes. We break the spherical symmetry of our 2D 
simulations by heating the initial profile with a combination of 75% spherical mode heating 
and 25% sin^ 6 mode heating. Thus for the 2D cases dQ/dtdV becomes 

dQ p f r - Tg 



dtdV 



(V)^^"^^ (1 + 0-5^^^'^) (2) 




10^° 



FIG. 1: The initial, t = 0, density profile used in our SN simulations. The dense core (A) inside 
20 km contains approximately 3Mq. The slow, outward-moving standing accretion shock (B) is 
located at 200 km. Above that we have a collapsing, 13.2 Mq, progenitor (C). 

In either case the total energy, Q, input during the simulation was recorded. Below the gain 
radius we wish to maintain the spherical density/gravity conditions and to prevent mass from 
escaping the core so that we do not disturb the shock heating and convective flow in the 
outer star. In a real supernova the core is stabilized by neutrino emission but in lieu of full 
implementation of this cooling mechanism our simulations achieve the same result by forcing 
all radial and angular velocities below Vg to zero. In this way we separate the evolution of 
the outer layers of the star from the complex behavior near the core. The runtime of the 
ID simulations is sufficiently short that we can vary the total energy deposition thereby 
obtaining a variety of results ranging from weak explosions through to the very powerful. 
The 2D simulations take much longer to run and so we have just one data set with an 
explosion energy set to 3 x 10^^ erg. 

Before presenting our simulation results we mention that in general, numerical schemes 
tend to spread the shock front over several zones and they therefore become artificially 



softened. Schirato and Fuller 27|] account for this artifact by steepening by hand the density 
profile at the shocks. The shocks in our simulations are sufficiently steep to demonstrate 
the primary effect on neutrino mixing, therefore we do not implement a similar correction. 
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B. Simulation Results 

The density profiles of the SN we obtain from the simulations possess noticeable differ- 
ences as the energy deposition changes. For weak explosions the profiles are monotonically 
decreasing functions of the radius, r, with a single, forward facing, shock front that moves 
relatively slowly outwards. As the deposition energy increases a lower-density cavity forms 
behind the shock. Further increases in Q eventually lead to the formation of a reverse shock 
behind the lower-density zone. All three features are present in the results of the 2D simu- 
lations but the higher dimensionality, coupled with the aspherical heating means, that this 
SN possesses a much more turbulent/chaotic profile. In figure ^ we show a time slice of 
the 2D SN simulation. The forward shock is seen at 'A', the 'reverse shock' is at 'B' and 
one of the many bubbles in the density profile is located at 'C. We note that although our 
applied neutrino heating is symmetric about the equator {0 = 7r/2) the resulting density 
profile is not. We attribute this to the inherent instability of multi-dimensional accretion 
shocks, potentially excited by simulation numerics. In what follows we discuss the behavior 
and properties of the various features of the profile. 

1. The Forward Shock 

The forward shock is a generic feature of supernovae simulations. A stalled forward shock 
was present in the initial profile inserted into the hydrodynamical code and the heating we 
introduced was meant to revive it's outward motion. After revival the shock propagates 
out through the star and is the explosion feature furthest from the proto-neutron star. The 
forward shock is visible in figure ([2]) as the large jump in density at larger radii and we 
note that the forward shock in the 2D simulation is slightly oblate due to the aspherical 
heating of the material. We study the behavior of the forward shock with various explosion 
strengths using the ID simulations. 

In Fig. ([3]) we show four snapshots of the density profile in a one- dimensional simulation 
in which the total energy deposition was Q = 1.66 x 10^^ erg. Actual SN are thought to 
be more energetic than this so we regard the results from this simulation as being at the 
lower end of realistic possibilities. This energy is also significantly less than that used in the 
simulation shown in Fig. ([2]) and, consequently, the profile is much simpler. The forward 
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FIG. 2: The density as a function of the radius and angle in a 2D SN model at t = 2.5 s. The 
forward shock is located to the left of 'A', the reverse shock is the step-up in density found to the 
right of 'B', and one of the many local cavities in the profile between the shocks is to the right of 
'C. 

shock is clearly visible and we note that the fractional jump in density across it does not 
vary. For normal shocks the density jump across the shock is given by 

Ap ^ (7 + m' 

p (7- l)M2 + 2' ^ ^ 
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FIG. 3: The density function of the radius in a ID SN model with Q = 1.66 x 10^^ erg at 
t = 0.9 s (dotted), t = 1.8 s (sohd), t = 3.6 s (long dashed) and t = 7.2 s (dash dot). 

where M the Mach number and 7 is determined by the ratio of specific heats. In our 
simulations we adopted a fixed equation of state and the Mach number of the shock is 
always sufficiently large that we find that the density jump does not vary to any great 
extent with time and/or with Q. But the position of the forward shock as a function of time 
is greatly influenced by Q with larger values leading to more rapidly moving forward shocks. 
This can be seen in Fig. (jlj) which are snapshots at the same moments as those in Fig. 
([3]) for a ID simulation with Q = 3.07 x 10^^ erg. This figure also displays the low-density 
cavity that can form behind the forward shock as Q increases. This profile is very similar 
to the profiles used by Fuller & Schirato [27I and Fogli et al. j2^. 




2. The Reverse Shock 

The heating that led to the regeneration of the forward shock continues to accelerate the 
material above the proto-neutron star even after the shock has been revived and is moving 
outwards. A wind is created with a velocity that increases with radius. When the velocity 
of the material becomes larger than the local sound speed a reverse shock is formed. The 
reverse shock feature was not present in the 'initial' profile: it develops only later. If a 
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FIG. 4: The density function of the radius in a ID SN model with Q = 3.07 x 10^^ erg at 
t = 0.9 s (dotted), t = 1.8 s (sohd), t = 3.6 s (long dashed) and t = 7.2 s (dash dot). 

reverse shock is to form enough energy must be deposited to create a sufficiently strong 
wind. A reverse shock can be seen in figure ([2]) where it is the step-up in density at the 
back of the turbulent zone behind the forward shock. The reverse shock is oblate due to the 
aspherical heating of the material in the simulation. 

Once again we can use ID simulations to study this feature more easily. No reverse shock 
was seen in the proffies shown figure ([3]) or (jlj) but when we increase Q to Q = 3.36 x 10^^ erg 
we obtain the results presented in figure ([5]). We notice that the forward shock is, again, 
moving more rapidly compared to the results shown in figures ([3]) and (jl]) due to the larger 
energy deposition Q. The figure also indicates that the reverse shock is smaller than the 
forward shock. As Q increases further the density jump across the reverse shock increases 
and it moves closer to the forward shock. We find that it is also possible for the reverse 
shock to penetrate to densities lower than the forward shock. Both of these behaviors 
are shown in figure ([6]) where we show snapshots of the profile for a ID simulation with 
Q = 4.51 X 10^^ erg. The density jump across the reverse shock is still smaller than across 
the forward but they are almost equivalent, and the density immediately behind the reverse 
shock is lower than the density immediately in front of the forward shock. It was mentioned 
by Tomas et al. [sO^ that this feature was also seen at early times in their simulations. But 
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r [cm] 

FIG. 5: The density function of the radius in a ID SN model where Q = 3.36 x 10^^ erg at 
t = 1 s (dotted), t = 1.5 s (sohd), t = 2 s (long dashed), t = 2.5 s (dash-dot) and t = 3 s (dash 
double dot). 




FIG. 6: The density as a function of the radius in a ID SN model with Q = 4.51 x 10^^ erg at 
t = 1.0 s (dashed), t = 1.5 s (solid) and t = 2.0 s (long dashed). 



12 



realistically SN may struggle to achieve explosions with comparable energies so we regard 
the results of this simulation as being at the upper end of possibilities. 

In both figures (jS]) and (jS]) the reverse shock was driven radially outward by a wind 
generated by the material heated above the proto-neutron star. Our energy deposition 
decreased exponentially with time which led to a gradual decrease in wind strength. In 
turn, as the wind abates, the outward motion of the reverse shock slows and, eventually, 
we find that its motion can be completely halted. This stalling of the reverse shock also 
occurred in the two-dimensional simulation with the shock stalling at different times for 
different polar angles. We also find that in both ID and 2D simulations the density jump 
across the reverse shock decreases as the shock is about to turn around. After the reverse 
shock stalls both the ID and 2D simulations indicate that the reverse shock feature then 
moves back towards the core. In figure ([5]) we see this stalling of the reverse shock and 
the backwards motion can be seen after comparing the profiles at t = 2.5 s and t = 3 s. 
The reverse shock in this simulation actually reached it's furthest radial position at about 
t = 2.4 s. Similar backwards motion for the reverse shock may be seen in figure (1) of Tomas 
et al. 3ol |. 

Eventually the reverse shock reaches the core whereupon our simulations indicate that 
it is reflected and subsequently becomes a weak forward shock. Further snapshots of the 
density profile taken from the simulation with Q = 3.36 x 10^^ erg are shown in figure 
where we see more clearly the backwards motion of the reverse shock and its later 
reflection. The radial position at which the reverse shock stalls depends upon the energy 
deposition. Additional snapshots from the simulation with slightly larger energy deposition, 
Q = 4.51 X 10^^ erg, are show in figure ([8]). Though the reverse shock in this simulation 
also attained it's maximum radial position at t = 2.4 s the figure shows that it was located 
further out before it was turned around. 

The reverse shock is an interesting feature of the SN and, in contrast with the forward 
shock, can move both outwards and inwards. For very energetic explosion it may penetrate 
to lower densities than that immediately in front of the forward shock. The density jump 
across the reverse shock can vary with time and does so most noticeably when the shock is 
about to stall. 
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FIG. 7: The density as a function of the radius in the same ID SN model shown in figure ([5|) i.e. 
with Q = 3.36 x 10^^ erg, at t = 2.8 s (dotted), t = 3.3 s (solid) and t = 4.3 s (long dashed). 
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FIG. 8: The density as a function of the radius in the same ID SN model shown in figure ([6]), that 
is (5 = 4.51 X 10^^ erg, at f = 4.0 s (dotted), t = 4.5 s (solid) and t = 4.9 s (long dashed). 
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FIG. 9: The density at a polar angle of 25° as a function of the radius in a 2D SN model at t = 3.9 s 
(solid), t = 4.8 s (dotted) and t = 5.7 s (dashed). 

3. Asphericity 

The difference in the hydrodynamics between the one and two dimensional SN simulations 
can be traced back to the aspherical heating of the material above the proto-neutron star 
in the two dimensional case. Non-spherical heating can lead to turbulent fluid flow creating 
eddies and bubbles as shown in figure ([2]). The neutrinos released by the proto-neutron 
star propagate along radial slices of the profile. In Fig. we show the density profile from 
the 2-D model at the polar angle of = 25° for snapshots at t = 3.9 s, t = 4.8 s and t = 5.7 s 
while in figure (fTOl) we present the density along three different polar angles at t = 2.5 s, 
i.e. taken from figure ([2]). In both we see the forward and reverse shocks and, as mentioned 
earlier, figure (fTOl) also shows that the radial position of both the forward and reverse shocks 
varies with the polar angle. Figure ([9]) also indicates that between t = 4.8 s and t = 5.7 s 
the reverse shock (along this radial slice) stalled and began to move back to the core. For 
other radial slices the reverse shock turnaround time will be different: in figure (ITT]) we plot 
radial slices along the 125° line of sight at various times. The reverse shock in the figure, 
located just beyond r = 10^ cm, turns around between 5.4 s and 6.3 s. 

In all our 2D figures the general shape of the density profile between the two shocks is 
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FIG. 10: The density as a function of the radius in a 2D SN model at t = 2.5 s. The angular slices 
are at 45° (solid), 105° (dotted) and 165° (dash dot). 
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FIG. 11: The density at a polar angle of 125° as a function of the radius in a 2D SN model at 
t = 4.5 s (dash dot), t = 5.4 s (solid), t = 6.3 s (dotted), and t = 7.2 s (dashed). 
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reminiscent of that shown in figure ([5]) for the ID simulation with a similar value for Q - 
i.e. a lower density region in front of the reverse shock, a wall of higher density material 
behind the forward shock - but clearly there are large 'fiuctuations' upon this basic trend 
both between the shocks and behind the reverse shock. 

C. Summary 

The forward shock, the reverse shocks, contact discontinuities and the local bub- 
bles/cavities are the features in the density profile with the greatest potential to alter the 
state of any neutrinos propagating through the SN since they represent the locations within 
the SN where the density gradient is largest. Large density gradients lead to non-adiabatic 
evolution of the neutrinos and significant differences compared to the neutrino propagation 
through the undisturbed profile. 

III. NEUTRINO MIXING 

The vast majority of the neutrinos emitted by the supernova last interact with matter 
at the neutrinosphere located at the surface of the proto-neutron star. As they propagate 
outwards a small percentage will be absorbed by nucleons thereby transferring energy to 
the supernova and, it is thought, reviving the stalled shock. But even though the remainder 
of the neutrinos survive the passage through the material overlying the proto-neutron star 
what emerges is not the same as what was emitted. The change in the neutrinos is due to 
neutrino oscillations and the presence of matter modulates this mixing. 

Neutrino oscillations arise due to a distinction between the interaction eigenstates (oth- 
erwise known as the fiavor eigenstates e, /x, r) and the eigenstates of the free Hamiltonian 
(known as the mass eigenstates) with masses mi, 7712,^3. Since there are three basis states 
a general neutrino wavefunction is described by three complex coefficients and evolves ac- 
cording to the Schrodinger equation. In the vacuum the Hamiltonian is diagonal in the 
mass basis but possesses off-diagonal terms in the flavor basis that are the cause of flavor 
oscillations. In the presence of matter a potential, V^(r), that takes into account coherent 
forward scattering of the neutrinos, is included in the Hamiltonian. Since we are only con- 
cerned with mixing between active neutrino flavors (i.e. all the flavors that have ordinary 
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weak interactions) we may subtract off tfie common neutral current contribution (wfiidi will 
contribute only an overall phase) leaving just the charged current contribution to the Ve — 
component of V{y). This contribution is the well-known Vee(r) = -\/2 G^Felr) p(r)/'"T'Ar 
where Gp is Fermi's constant, p(r) is the mass density, Yf.{Y) is the electron fraction and 
mjv is the nucleon mass. The effect of the matter upon the antineutrinos differs from that 
of the neutrinos - the potential V{v) that appears in the antineutrino Hamiltonian has the 
same magnitude but the opposite sign. Due to the inclusion of T4e(r) the Hamiltonian is a 
function of position and is neither diagonal in the mass basis nor the flavor basis. One may 
try and diagonalize H{r) but the unitary transformation that relates the flavor basis to the 
new 'matter' basis also varies with the position. Consequently the gradient of the unitary 
transformation is non-zero and one finds that the Schrodingcr equation in this new matter 
basis - that was meant to diagonalize H - picks up off-diagonal terms. Thus it is not possible, 
in general, to find a basis that diagonalizes the Hamiltonian and therefore oscillations of the 
amplitudes that describe the wavefunction occur in every basis. 

In general, the three complex components of the wavefunction oscillate simultaneously. 
Exactly what occurs depends upon the energy of the neutrino the differences between the 
squares of the masses 1711,1712,1713, the mixing angles that describe the relationship between 
the flavor and mass basis, and, of course, Vee(r). Since there are three mass eigenstates there 
are three separate mass splittings Sm^j — ml — (though only two are independent) and 
the relationship between the neutrino flavor and mass bases is described by the matrix U 
which, in turn, is parameterized by three mixing angles ^12, ^13 and ^23 plus a CP-phase 5. 
The structure of U is 



U 



^ C12 Ci3 S12C13 Si3e ^ 

"•512 C23 ~ C12 Si3 S23 e*'' C12 C23 — S12 Si3 S23 e"^ Ci3 S23 



(4) 



y 512 S23 - C12 Si3 C23 e'" -C12 S23 - Si2 Sl3 C23 e-" Ci3 C23 J 

where Cij = cos 6 ij, Sij = sin 9ij. Mixing has been observed in the neutrinos emitted by 
the Sun and the neutrinos produced by cosmic rays striking the atmosphere. Both have 
been confirmed with terrestrial experiments. Each observation of neutrino mixing can be 
described by a single Sm^ - pair of parameters and experimentally the 'solar' and 'at- 
mospheric' mass splittings differ by around a factor of ~ 30 with the solar mass splitting, 
5mQ, being the smaller of the two. This observation permits us to consider the evolution 
of the general, three component, neutrino wavefunction as being factored into spatially dis- 
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tinct, two- neutrino mixes. The factorization simplifies matters greatly and from it one can 
demonstrate that there are two resonances in the supernova density profile: the so-called 'L' 
resonance and the 'H' resonance. The L resonance, at lower density and thus further from 
the proto-neutron star, involves mixing between matter states ui and U2, the relevant mass 
splitting is Sm^i - which is approximately the solar mass splitting 6mQ - and the mixing 
angle is 612 - which is approximately the mixing angle determined by the solar neutrino 
experiments 6q. The relevant mass eigenstates and mass splitting for the H resonance, at 
higher density and closer to the core, depend upon the sign of (5m|2 (or, equivalently, Smli) 
and this is not currently known. If 6771^2 is positive (a normal hierarchy) then the H res- 
onance involves mixing between states 1^2 and 1/3 and the mass splitting If f^^L 
negative (an inverted hierarchy) then at the H resonance it is the antineutrinos states ui and 
z/3 that mix and the relevant mass splitting is 5m|]^ = 577132 + 5^21. In either case the mixing 
angle is the unknown, but small, ^13. Quite generally the coherent matter basis wavefunction 
that arrives at a distance, d, from the core of the SN is related to the initial matter basis 
wavefunction at the proto-neutron star via the equation ipuid, E) = Sy{d, E) V^i,(0, E) while 
for the anti-matter states ipy{d^E) = Sp{d, E) ipp{0, E). The matrices Si, and Sp are the 
S-matrices for the neutrinos and antineutrinos respectively, for a discussion of this approach 
to neutrino oscillations, see e.g. j2o|. To determine ip^{d,E) and iljp{d,E) we need to know 
both the initial states ipu{0, E) and ^/'p(0, E) and the two matrices and Sp. 

The matrices S^{d,E) and Sp{d,E) can be factored as S^{d,E) = Sy^ Sl,u{E) Sh,u{E) 
and Sp{d,E) = Sv,u Sh,u{E) which are sufficiently general to accommodate our lack of 
knowledge of the hierarchy. The two matrices Sv,u and Sv,p represent the neutrino propaga- 
tion from the surface of the SN through the vacuum to Earth and both are diagonal in the 
mass basis. The matrices Sl,u{E), Shm{E) and Shm{E) represent the change to the initial 
wavefunction due to the neutrino's or antineutrino's passage through the H or L resonance^. 
There is no matrix Sl^p{E) because antineutrinos do not experience an L resonance in mat- 
ter. The two matrices Sh,v{E) and Sh,d{E) are evaluated somewhere between the H and 
L resonances while Sl{E) is evaluated at the surface of the SN. These evaluation positions 

^ Many of the profiles presented in section S|TT]havc multiple H resonances therefore the matrices Sh,u{E) 
and Sh.v{E) represent their combined effect of all the H resonances and Sl^^{E) all the L resonances. We 
assume that the neutrinos encounter all L-type resonances after all H-type resonances. This is the case 
for all profiles we have generated 
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will not affect our result if they are sufficiently far from the resonances. The structures of 
Sl,v{E), Sh,v{E) and Sh,d{E) follow from our knowledge of the matter states that mix at 
either the H or L resonance and are thus 
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after omitting irrelevant phases. Prom these matrices we define 

PAE) = 1-ME)\'^\ME)\', 
Ph{E) = I - \aH{E)\^ = \h{E)\\ 
PAE) = l-\aAE)?^\ME)\\ 



(8) 
(9) 
(10) 



which are the crossing probabilities for neutrinos or antineutrinos at the two resonances. 
The resonances are said to be 'adiabatic' or 'non-adiabatic' depending upon whether the 
crossing probability is close to zero or closer to unity. Indeed these are the two natural 
values since the crossing probability is determined by the ratio of the resonance width to 
the local oscillation lengthscale and typically one is significantly larger than the other. 

Our interest now turns to the initial states. The density at the proto-neutron star 
is so large that the matter cigenstates and the flavor eigenstates are strongly aligned 
there. From a full 3-neutrino mixing calculation we find that the initial matter basis 
spectra for a normal hierarchy (NH) are $^,3(0, £■) = ^„^{Q,E), $p^(0,£') = $,7^(0, £■) 
while all other states, ^ui{^,E),(^A^^E),(^i^^{Q,E) and $,73(0, £■) are equal to the 
$^^(0,£;),$^,(0,^),$p^(0,^),$p,(0,E) spectrum which we call $^,(0,E). With an in- 
verted hierarchy (IH) the initial states arc $1,2(0, -E) = $zy^(0,i?), ^p-A'^^E) = $£7^(0, E) and 
this time ^^.^(O, E), $^3(0, E), $^,,(0, E) and ^^.^(O, E) are all equal to $^^(0, E). 
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Putting the initial spectra and definitions of tlie S-matrices togetlier one finds that the 
flux of matter state i a distance d from the supernova is given by 



F^ME) = ^ ^|(5.),^.(rf,i?)r$.^,(0,i?), (11) 

j 

F.ME) = J2\(^^h(d^^)\'^^ME)- (12) 

j 

where (0,-E) and $pj(0,i?) are the initial spectra of the matter states. But what are 
detected at Earth are, of course, the flavor states. During their flight from the supernova 
to Earth any coherence between the matter eigenstates is lost so the flux of flavor a that 
arrives at Earth is the incoherent sum 

F,ME) = J2\Uai\'F,Xd,E). (13) 

i 

These flavor fluxes may be rewritten in terms of the fluxes emitted by the neutrinosphere 
by introducing the survival probabilities p{E) and p{E) for the electron neutrinos and an- 
tineutrinos respectively since these are the two flavors with distinct initial spectra. Written 
this way the detectable flavor fluxes are 

F,ME) = j^ [piE)^,ME) + il-pm^uME)], (14) 

F,ME) = j^ [piE)^,ME) + il-piE))^,ME)], (15) 

AF^M E) = [(1 - pm $,^(0, E) + {1- p{E)) $p^(0, E) + {2 + p{E) + p{E)) $.,(0, E)] 

(16) 



where we have adopted the notation of Dighe & Smirnov [13(] by denoting F^j^ +F^^+Fp^+Fp^ 
by 4Fj,^. From equations f|TT]) and f|T2|) . and the initial spectra, we flnd that p{E) and p{E) 
are related to the elements of S{E) and Sp{E) via 

PiE) = J2\Ue^\'\{Su\,{E)\' (17) 

i 

p{E) = Y.\Ue^\'\iS,UE)\' (18) 
i 

where j = 3, k = 1 for a normal hierarchy, and j = 2, k = 3 for an inverted hierarchy. After 
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inserting the definition of tlie crossing probabilities p{E) and p{E) become 



NH 



IH 



p{E) = |f/eip Pl{E) Ph{E) + |f/e2p (1 " Pl{E)) Ph{E) + |f/,3p (1 " Ph{E)) 

p{E) = |f/eip (1 - Ph{E)) + |f/e3p Ph{E) 

(19) 

p{E) = \Uel\^ Pl{E) (1 - Ph{E)) + |f/e2p (1 " Pl{E)) (1 - Ph{E)) + |f/e3p Ph{E) 
p{E) = It/elP Ph{E) + |f/e3p (1 - Ph{E)) 

(20) 

Tliese formulae are quite general but from our knowledge of the various neutrino oscillation 
parameters we can be more specific about what exactly happens to the neutrinos and an- 
tineutrinos. The solar mixing parameters are such that the neutrino L resonance is almost 
completely adiabatic until, perhaps, the very end of the supernova neutrino signal. For this 
reason Pl{E) = 0. Thus at the end we find that determining the neutrino fluxes at Earth 
comes down to computing Ph{E) or Ph{E). 

The Schrodinger equation forms a starting point by which Ph{E) or Ph{E) can be 
determined. But if one naively applies, for example, a Runge-Kutta integrator to this 
equation one quickly discovers that this is a difficult problem from a numerical standpoint 
because such algorithms are simply not suited to differential equations where the solution 
is a highly oscillatory function. A number of alternate methods have been developed for 
calculating Ph{E) or Ph{E), such as application of the the Landau-Zener result or the 
semi-analytic method by Balantekin & Beacom [35], but for one reason or another these 



alternate approaches can break down for complex profiles. Some authors 20|, l36|, |37|, |38| have 



recognized that the evolution of the neutrino wave can be recast as a scattering problem and 



we adopt in this paper the algorithm outlined in Kneller & McLaughlin 20|] which computes 
the evolution of the neutrino wavefunction by a Monte Carlo integration. For this paper we 
selected \6m^\ = 3 x 10~^eV^ and for ^13 we use sin^ ^^13 = 10~^ as a representative value. 
With this choice of 6*13 the unperturbed profile is just sufficiently adiabatic to suppress 
the transformation effects of all but the most significant density features. The current 
experimental limit on ^13 is set by the CHOOZ experiment jsQ], and for our choice of 
that limit is sin^ ^13 < 0.1. 

Finally, in addition to the effects caused by the propagation through the supernova there 
are potentially Earth matter effects that can arise. These are straightforward to compute 
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and we do not consider them here since their presence (or absence) in the neutrino signal 
is a function of the position of the supernova with respect to the detector when the event 
occurs. 

IV. PROFILE FEATURES AND THE EFFECTS UPON THE CROSSING PROB- 
ABILITY 

With the SN simulations complete and a method for calculating the crossing probabilities 
in hand we turn to study the effects of various features in the profiles upon Ph both as a 
function of time and energy. 

A. The Forward Shock 

The forward shock is a generic feature of all our SN simulations. Initially the shock is 
located close to the core and at high density and then, after it is regenerated, begins to 
move outwards and to lower densities as shown in the figures from section ^|TT1 In Fig. f fT2|) 
we reproduce the results from Fig. (JSj) and also superimpose the resonance densities for 5, 
10, 20, 40 and 80 MeV neutrinos. The radii where the horizontal dashed lines intercept 
the profile are the locations of the resonances. And in Fig. f|T3l) we show Ph as a function 
of neutrino energy at six snapshots of this same simulation. The crossing probability as a 
function of the energy clearly evolves with time. If ^13 is not too small then the evolution 
of the neutrinos through the undisturbed progenitor profile is almost adiabatic i.e. Ph ~ 0. 
As the forward shock arrives at the H-resonance for 5 MeV neutrinos - Fig. (fT2i) indicates 
this occurs at shortly after t ~ 1.8 s - the evolution becomes non-adiabatic i.e. Ph ~ 1 
because the density jump across the shock straddles the resonance densities of this neutrino 
energy. As time progresses and the shock moves outwards to lower densities we see from 
Fig. (fT2|) that the shock will begin to affect 40 MeV neutrinos at t = 3.6 s. The non- 
adiabaticity sweeps up through the neutrino spectrum from low energy to high. As time 
progress further eventually the shock ceases to affect the neutrinos of a particular energy 
and so their propagation returns to being adiabatic. For this simulation Fig. ([3]) indicates 
that the evolution of 5 MeV neutrinos returns to being adiabatic at around t = 3.6 s and 
that for 40 MeV the return occurs roughly at t = 7.2 s. Note that the 5 MeV neutrinos 
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FIG. 12: The density as a function of the radius in a ID SN model with Q = 1.66 x 10^^ erg 
at t = 0.9 s (dotted), t = 1.8 s (sohd), t = 3.6 s (long dashed) and t = 7.2 s (dash dot). The 
horizontal dashed lines are (from top to bottom) the resonance densities for 5, 10, 20, 40 and 
80 MeV neutrinos. 

were affected for a much briefer period (~ 1.8 s) than the 40 MeV neutrinos (~ 3.6 s). This 
is due to the shape of the density profile. A complimentary perspective is to look at how 
particular neutrino energies evolve with time such as those in Fig. ( JT^ . The temporary 
transition to non-adiabaticity for each neutrino energy is clearly visible in the figure. This 
figure makes it most obvious that the lower energies are affected before the higher and also 
that the duration of the non-adiabatic period increases with the neutrino energy. 

The extent of the shock feature in neutrino energy seen in Fig. f[T^ is related to the 
density jump, Ap, across the shock. At any given time, if Es is the highest neutrino energy 
affected by the shock and AE is the range of neutrino energies then the density jump across 
the shock is 

An Af" 

(21) 



Ap 
P 



- AE 



which is independent of the mixing parameters. From equation ((31) we saw that, for normal 
shocks, Ap is a function of the Mach number and ratio of specific heats 7. If M is large 
then Ap, and consequently AE, are essentially just a function of 7 but if the Mach number 
M is not too large - which may be the case - then the width of this feature could be used to 
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FIG. 13: The H resonance crossing probability Ph as a function of neutrino energy for the ID SN 
simulation with Q = 1.66 x 10^^ erg. In the top panel the snapshot times are t = 2 s (solid) and 
t = 2.5 s (dashed), in the middle panel t = 3 s (solid) and t = 3.5 s (dashed) and in the bottom 
panel t = 4 s (solid) and t = 4.5 s (dashed). 

infer M if 7 is known. As the SN proceeds to explode Es moves up through the neutrino 
spectrum. The correspondence between energy and resonance density means that if we 
know, or assume, a progenitor profile shape then Es measures the shock position rs- From 
the locations of Es shown in Fig. (fT3|) . and the initial profile, we show in Fig. (ITSl) rs 
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FIG. 14: The H resonance crossing probability Ph as a function of time for selected neutrino 
energies through the ID SN model with Q = 1.66 x 10^^ erg. The curves are: E = 10 MeV 
(solid), E = 15 MeV (long dashed), E = 20 MeV (dash-dot), E = 25 MeV (short dashed), and 
= 30 MeV (dash double-dot). 

as function of time for this particular simulation. The change of with time of course 
measures the shock speed but we can also use rsit) to extrapolate the shock position both 
forward in time, if necessary, but, more interestingly, also backwards toward the proto- 
neutron star. If this extrapolation can be accomplished successfully then we can infer the 
moment, ^200, when rs ~ 200 km. The neutronization burst (if detected) supplies a zero of 
time so t2oo 7^ is the brief period the shock was stalled. Thus the neutrino signal, and in 
particular the location of the forward shock, contains evidence of whether a key component 
of the core-collapse supernova paradigm is correct and can provide quantitative data with 
which to compare with more sophisticated SN simulations than ours. We illustrate this idea 
in Fig. (fT5|) where we have artificially included an offset in time. 

B. The Reverse Shock 

The reverse shock was a feature seen in the more powerful explosions from section ^|TT] 
and the 2D result. This feature was not present in the initial profile and appeared later 
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FIG. 15: The shock position as a function of time, rs{t), for the ID SN simulation where Q = 
1.66 X 10^^ erg. An artificial time delay of 500 ms has been added to mimic the stalling of the 
shock at rs ~ 200 km. 

on when when the velocity of the wind superseded the local sound speed. After its forma- 
tion the reverse shock initially moved outwards into the star but, as the energy deposition 
faded with time and the strength of the wind abated, eventually the reverse shock stalled 
and then headed back to the core. Like the forward shock, the reverse shock affects the 
adiabaticity of the evolution of the neutrino wavefunction through the profile so that neutri- 
nos with resonance densities straddled by the density jump across the reverse shock evolve 
non-adiabatically as they pass through it. Note that any neutrino with a resonance density 
straddled by the density jump across the reverse shock will also experience two other res- 
onances: one before the reverse shock and one after. This is shown in Fig. (fT6|l where we 
have superimposed the resonance densities for 5, 20, and 80 MeV neutrinos upon the profiles 
shown in Fig. (Q. We see from the figure that it is possible for some neutrinos to be affected 
by both shocks. If this occurs, and if we neglect the effects from other resonances in the 
profile, then we might expect that the net effect upon such doubly shock-affected neutrinos 
to be zero 30|. However this expected cancellation is not seen in our results for Ph as a 
function of the neutrino energy shown in Fig. f|T7j) for the model where Q = 3.36 x 10^^ erg. 



The crossing probabilities plotted in the figure show some similarities to those plotted in Fig. 
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FIG. 16: The density as a function of the radius in a ID SN model with Q = 4.51 x 10^^ erg at 
t = 1.0 s (dashed), t = 1.5 s (sohd) and t = 2.0 s (long dashed). The horizontal dashed lines are 
(from top to bottom) the resonance densities for 5, 20 and 80 MeV neutrinos. 

([T3|) at the higher neutrino energies affected by the forward shock where there is change from 
adiabatic to non-adiabatic evolution. But for those energies where cancellation is naively 
expected we see instead that Ph oscillates wildly. These rapid oscillations are phase effects 
due to the interference between the two shocks. Similar rapid oscillations in the crossing 



probability have been seen previously in SN neutrino calculations by Fogli et al. [29|] in a 
profile with a forward shock and then a bubble cavity behind it, in the SN test case used 
by Kneller & McLaughlin |2o|, and in the results of Dasgupta & Dighe 21 1 where their 
significance was emphasized. 

Fig. ( JTSil we display the crossing probability for the model where Q = 4.51 x 10^^ erg at 
t = 1.1 s, t = 1.4 s and t = 3.0 s. Again phase effects are seen. Compared to the results 
shown in Fig. (11 7p for the weaker explosion where Q = 3.36 x 10^^ erg and in Fig. (ITB]) 
for Q = 1.66 X 10^^ erg, in this case there is no indication of the characteristic transition 
from adiabatic to non-adiabatic propagation associated with the forward shock. This occurs 
because, as noted earlier about this model, the reverse shock penetrates to lower densities 
than the forward shock at these early times. For this particular model the reverse shock 
affects the adiabaticity of a particular neutrino energy before the forward shock. 
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FIG. 17: The H resonance crossing probability Ph as a function of neutrino energy for three 
snapshots taken from the ID SN model where Q = 3.36 x 10^^ erg. In the top panel the time is 
t = 1.8 s, in the middle t = 2.0 s, and in the bottom panel t = 2.4 s. 

The presence of phase effects in Ph are caused by interference between forward and 
reverse shocks. But it would be nice to find in the neutrino signal a clean signature of the 
reverse shock that is not contaminated by the forward shock or other features of the profile. 
A potential signature could arise from the behavior noted earlier which is that the reverse 
shock stalled as the energy deposition into the material above the proto-neutron star faded 
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FIG. 18: The H resonance crossing probability Ph as a function of neutrino energy for three 
snapshots taken from the ID SN model where Q = 4.51 x 10^^ erg. In the top panel the time is 
t = 1.1 s, in the middle t = 1.4 s, and in the bottom panel t = 3.0 s. 

and then headed back towards the core. As the reverse shock moves backwards its effects 
move down through the neutrino spectrum and might eventually cease to overlap with the 
forward shock. Since the density jump across the reverse shock becomes quite small by the 
time this behavior occurs the reverse shock will produce a narrow spectral feature. For the 
simulation where Q = 3.36 x 10^^ erg, shown in Fig. Q, this expected pattern is difficult 
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to see in the neutrino signal because the forward shock has not swept through the neutrino 
spectrum by the time the reverse shock turns around. But in the simulation with slightly 
larger energy deposition, Q = 4.51 x 10^^ erg, the forward shock moves much more quickly 
and has largely swept through the H resonance region by the time the reverse shock starts 
to make its way back to the core. For this simulation the backwards moving reverse shock 
is more visible in the neutrino signal. The crossing probability for this simulation at the 
snapshots shown in Fig. ([8]) are presented in Fig. f|T9l) . The backwards moving reverse shock 
is responsible for the phase effects at E^, ~ 10 MeV in the middle panel. And then once the 
reverse shock has reached the core our boundary conditions there led to its reflection and 
so the reverse shock became a outward moving, weak forward shock and re-appeared in the 
neutrino spectrum as a spectral feature moving from low to high neutrino energies. This is 
the cause of the bump in Ph seen at ^ 15 — 20 MeV in the bottom panel. 

In summary, though it is apparent that the neutrino signal can vary considerably de- 
pending upon exactly how the reverse shock behaves the presence of strong phase effects - 
rapid oscillations with large amplitude in both time and energy of the crossing probability 
- are a notable signature of the presence of multiple shocks. Therefore, there exists tremen- 
dous potential in a future detection of the a supernova neutrino signal for probing the inner 
hydrodynamics of the explosion. 

C. Asphericity 

The notable difference between the one-dimensional and the two-dimensional models was 
the local density fluctuations that appeared in the latter. These local density enhance- 
ments/cavities affect only neutrinos along particular lines of sight so if one could determine 
their presence in the signal we would obtain an indication of the asphericity of the explosion. 
The most obvious method would be to compare different lines of sight to the same SN but 
obviously this is not practical. Instead we must find another signature of asphericity that 
one might use. In figure (!20|) we show the crossing probability Ph as a function of neutrino 
energy for the snapshot t = 2.4 s, t = 5.4 s and t = 6.4 s and a line of sight 6 = 25°. 
Features in the figure resemble those seen in the ID results: the forward shock leads to a 
range of energies with Ph ^ 1 as seen in the top panel, phase effects due to presence of 
the reverse shock lead to the high 'frequency' oscillations of Ph seen in the middle panel. 
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FIG. 19: The H resonance crossing probability Ph as a function of neutrino energy for the ID SN 
model where Q = 4.51 x 10^^ erg. From top to bottom the snapshot times are t = A s, t = 4.5 s 
and t = 4.9 s. 

and the turn around of the reverse shock can be seen in the bottom panel at ~ 15 MeV 
where it appears, as before, as a narrow range of high frequency phase effects moving down 
through the spectrum. While a detailed analysis of these results might indicate a difference 
that is due to asphericity there is no striking feature that one can point to. However, at 
late times we found a change in the qualitative appearance of Ph versus neutrino energy. 
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FIG. 20: The H resonance crossing probability Ph as a function of neutrino energy for a radial 
slice at 6* = 25° through the 2D SN model. From top to bottom the snapshot times are t = 2.4 s, 
t = 5.4 s and t = 6.4 s. 

The high frequency phase effects no longer dominate the entire spectrum and now lower 
frequency portions show up. We did not see anything comparable at late times in any of 
our ID simulation results. The changing quality of the phase effects are seen in Fig. f l2T|) 
which shows the crossing probability at the H resonance as a function of neutrino energy 
for the same line of sight at = 25° through the 2D model. High frequencies can be seen 
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FIG. 21: The H resonance crossing probability Ph as a function of neutrino energy for a radial 
slice at 6* = 25° through the 2D SN model. From top to bottom the snapshot times are t = 7.4 s, 
t = 8.0 s and t = 9.0 s. 

above ~ 45 MeV in the top panel and above ~ 55 MeV in the middle panel but they 
have disappeared hj t = 9.0 s. In all three panels we see that at lower neutrino energies the 
curve is smoother with only narrow patches of high frequency phase effects that moved up 
through the spectrum. The shocks have largely ceased affecting the lower neutrino energies 
by this time so it is the 'fluctuations' upon the profile - both between the shocks and behind 



34 



the reverse shock - which produce mildly adiabatic resonances, and/or resonances that are 
located very close to one another (and possibly even overlap) that give rise to these features 
in the figure. 

The effect upon neutrino flavor transport of supernovae density fluctuations was consid- 
ered by Loreti et al [40] and later by Fogli et al 4l| and Friedland & Gruzinov [42]. These 
authors suggest that Ph should tend to 1/2 in the presence of many small scale fluctuations. 
We note, however, that due to the exponential growth of the radial grid spacing the scale 
of features we can resolve by r ~ 10^ km — 10^ km is of order 5r ~ 100 km and this is 
somewhat larger than the scale considered by Fogli et al and Friedland & Gruzinov. In 
our formulation of neutrino mixing using S-matrices we can motivate this expected result 
of Ph 1/2 by dividing the density proflle into domains such that within each there is a 
single neutrino resonance. The passage of a neutrino (or antineutrino) through each interval 
is described by an S-matrix which has exactly the same structure as equations ([6]) or ([7]) as 
appropriate. The magnitudes of the elements a and /3 he between zero and unity depending 
upon the adiabaticity of the resonance in that region. The multiplication of all these S 
matrices together (in the correct order) to achieve the net effect of passing through all the 
resonances will tend to produce a net result where either an and (3h or an and (3h are of 
equal magnitude. In our simulations we had cases where we found up to 13 resonances but 
have not been able to reproduce the anticipated limiting behavior. Nevertheless, as shown in 
Figs. ( !20l) and (!2T1) . we see a clear distinctive qualitative change in the crossing probabilities 
at late time in aspherical models. 



V. DETECTOR SIGNALS 

Since the detection of neutrinos from SN 1987A in 1MB [51] and Kamiokande [52;] the 
number and size of the neutrino detectors have both increased in scale so that many thou- 
sands or even tens of thousands of neutrinos will be detected from the next Galactic super- 
nova. Though many neutrino detectors are water Cerenkov detectors that are best suited 
to detecting the electron antineutrinos from supernovae via the inverse-/^ reaction on pro- 
tons, other neutrino detection technologies have been demonstrated or explored that allow 
one to see the electron neutrinos and the /i/r flavors via nucleus- neutrino reactions such as 
0{i'e,e)F, or neutral current reactions such as D{i'x,i^x)np. The detection of these other 
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neutrino flavors would provide a great deal of valuable information we could then use to 
learn about the SN explosion. For the purposes of demonstrating the temporal evolution 
of supernovae neutrino signals we shall consider two detector types: the water Cerenkov 
detector because the majority of neutrino detectors are of this type, and a heavy water 
detector (such as was used in the Sudbury Neutrino Observatory) since this detector was 



also considered by Schirato & Fuller 27|] . The largest neutrino detectors are water Cerenkov 



detectors. These detectors are sensitive to both electron neutrinos and electron anti neu- 



trinos via the reactions p(z/e, e)n and 0(z/e, e)F [43|| with the inverse-/? reaction dominating 
the signal. The cross sections for the neutrino-nucleon reactions are well known but a good 
deal of uncertainty exists in the important neutrino-nucleus cross sections. A heavy water 
detector is sensitive to both charged-current and neutral current events so could be used 
to detect the fi and r flavors. However obtaining energy resolution for the neutral current 



events is difficult so the observable quantity considered by Schirato & Fuller j27j was the 
ratio of total charged-current to neutral current events. Though temporal evolution of this 
ratio is somewhat smothered due to the fact that the event rates are integrated over the 
neutrino spectrum Schirato & Fuller demonstrated that variations of the ratio would occur 
at the level of 10 — 20%. Lead detectors are also capable of detecting the fi/r flux and have 



been studied by Fuller, Haxton & McLaughlin [12] Engel, McLaughlin & Volpe [14| and 
Kolbe et al £|. 

To determine what any given detector will see we must fold in three effects: the initial 
neutrino spectra emitted by the proto-neutron star, the energy dependence of the cross 
sections in the detector and, thirdly, the energy resolution of the detector. 

A. Initial Spectra 

There is a considerable degree of uncertainty about the emitted neutrino spectra. Fur- 



45| . Pastor, Raffelt & Semikoz 



^hermore, neutrino self interactions - described in Samuel 

46 ] and Friedland & Lunardini [l^ - in the region immediately above the proto-neutron 



star can alter the spectra after it's emission from the proto-neutron star but before the flux 
reaches the H resonance layer. There is a possibility that in some cases the spectra may be 
'equilibrated' by the time they get to the H resonance jl?! . In such circumstances the 
changes in the density proflle and the neutrino mixing described in previous sections would 
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not lead to any observable changes in a detector except for the gradual fading of the signal. 
We shall consider the case where such equilibration does not occur and experiment with 
an 'alpha fit' power law for the spectrum exiting the neutrinosphere motivated both by its 
analytic simplicity and by recent supernova core Monte Carlo simulations 3]. For the time 
dependence of the luminosity we use an exponential decay with a time constant of r = 3 s. 
The emitted differential spectra are then 



-(a + l) 



exp[-t/r]. (22) 



where L^, is the luminosity, (E^) the mean energy and a is the 'pinch' parameter. We adopt a 
hierarchy of mean energies: {E,^J = 12 MeV, {EpJ = 15 MeV and {E,y^) = 18 MeV (where, 
again, x stands for either fi,T,fL or f) and the luminosities are L^^ = Lp^ = 6 x 10^^ erg/s, 
L^^ = 2 X 10^^ erg/s. In Figs. ( l22l) and ( l23l) we present the neutrino spectra that would be 
observed at Earth from a SN at a distance of 10 kpc if the evolution were either completely 
adiabatic or completely non-adiabatic for example through a single H resonance, i.e. Ph = 
or Ph = 1 in either the inverted or normal hierarchy. The two curves in each panel of the 
figures therefore represent extremes in the sense that the actual flux at a given energy would 
fall somewhere in between the two curves. We do not plot the electron antineutrino flux for 
the normal hierarchy or the electron neutrino flux for the inverted hierarchy because these 
fluxes are always adiabatic and would only fade exponentially with time. In each panel of 
the two flgures the pinch parameter was set to a = 3, other than a change in shape the 
general behavior of the curves at other values of a is unchanged. There are many points 
to notice in each flgure. The flrst is that for any given neutrino flux and hierarchy there 
is a particular energy where the two curves cross. This will occur as long as the average 
energies of flavor neutrinos are greater than that of the Ue and z/g. The crossing can be 
seen best in the spectrum in the normal hierarchy. We deflne the energies where the 
curves cross as the critical energies. The critical energies vary with both flavor and the 
hierachy. At the critical energy for a particular hiearachy and flavor the adiabaticity of the 
neutrinos as they propagate through the proflle is irrelevant and the flux of that flavor and 
at that energy through a detector is independent of Ph- The other signiflcant feature to 
notice is the large difference between the two curves for the electron neutrinos and electron 
antineutrinos for energies below ~ 30 MeV. This wide separation of the extremal fluxes 
at low to medium neutrino energies, E < 30 MeV, suggests that we focus our attention 
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FIG. 22: Fluxes of electron antineutrinos (top panel) and the sum of i'^ + v'fj, + I'r + i^r neutrino 
fluxes (bottom panel) for completely adiabatic, i.e. Ph = (solid) and non-adiabatic, i.e. Ph = 1 
(dashed) evolution. The hierarchy is inverted, the pinch parameter was set to a = 3 and the 
vacuum mixing angles are sin^(2^i2) = 0.8, sin^(20i3) = 4 x 10~^. The distance to the SN was set 
at d = 10 kpc and the time t to t = 0.0 s. 
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FIG. 23: The fluxes of electron neutrinos (top panel) and the sum of i'^ + v'fj, + I'r + i^r neutrino 
fluxes (bottom panel) for completely adiabatic, i.e. Ph = (solid line) and non-adiabatic, i.e. 
Ph = 1 (dashed) evolution. Here the hierarchy is normal but otherwise all parameters are the 
same as those used in Fig. (f22|) . 

upon this portion of the spectrum because it is here that changes in adiabaticity will lead 
to the most significant changes in the flux. The flux at higher energies will also change with 
the evolving adiabaticity but to a lesser degree because the fluxes are so small. Thus if the 
evolution of the neutrinos evolves from being adiabatic to non-adiabatic then, in the case 
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of a normal hierarchy, the flux of electron neutrinos will increase for those energies below 
the critical energy and decrease above. Similar behavior will occur for antineutrinos in the 
inverted hierarchy. The decrement above the critical energy is difficult to see in the figures 
because the critical energies are all far into the tail of the spectra. For the + t'^ + z^r + 
neutrino fiuxes the opposite change occurs: a transition to non-adiabaticity leads to lower 
fiuxes of these fiavors at energies below the critical energy and higher fiuxes above. 

B. General Considerations 

Having determined the effects of the evolving density profile upon the neutrino mixing we 
are now in a position to examine the actual neutrino signal within a detector. As discussed 
in section llVt the actual neutrino fiux of a particular fiavor through a detector differs from 
that emitted from the neutrinosphere. In the simplest scenario, if the neutrinos propagate 
through the outer layers of the SN adiabatically then the survival probabilities p and p are 
found by setting Ph = Pl = Ph = in equations ( JT9l) and ( l20l) . For the normal hierarchy 
P = |f^e3p and p = |f/eiP- for the inverse hierarchy p = |t/e2p and p = |t/e3p- Using the 
value of sin^ 29i3 = 4 x 10~^ we used to calculate Ph and sin^ 29i2 = 0.8 based upon the 
solar mixing parameters, we see that for the normal hierarchy the amount of the emitted 
flux present in the electron neutrino flux at a detector here on Earth is minuscule. Similarly 
we see that for an inverted hierarchy the amount of the emitted electron antineutrino flux 
present in any detected electron antineutrino flux is equally small. Thus if the mass hierarchy 
is normal, sin^ 29is is not too small, and the neutrino propagation adiabatic then the electron 
neutrino detectors on Earth will essentially be detecting the emitted z/^ and Ur flux. Or if 
the hierarchy is inverted, again sin^ 2^13 is not too small, and the neutrino propagation 
adiabatic then the electron antineutrino flux at Earth is the emitted z/^ and Ur- In order to 
see the emitted electron neutrino or antineutrino flux we would need to detect the z/^, u^, 
Vj. and/or v^- flavors or seek the small amount of the emitted z/e/t'e present in the detected 
z/g/z/g flux. In either case, it will be very difficult to reconstruct the the z/g or z/g spectrum 
that was originally emitted from the proto-neutron star. 

The evolution of the density profile changes this conclusion. Depending upon the scenario, 
either p (for a normal hierarchy) or p (for an inverted hierarchy) will change with time and 
so the proportions of the emitted fiuxes present in the fiuxes at Earth will also change. 
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Note that p and p do not both evolve with time within in the same scenario; one will stay 
constant. This means there is a potential to separate the time evolution of the emitted flux 
from the neutrinosphere from the evolution of the density profile through which the neutrinos 
propagate. As the star explodes and the forward shock eventually reaches the H resonances 
the passage of portions of the neutrino spectrum through the density profile becomes non- 
adiabatic. Whether it is the neutrinos or antineutrinos that are affected is determined by 
the hierarchy. When the passage of the neutrinos is completely non-adiabatic we see that p 
or p change to p = |f/e2p or p = |f/eiP- Again using sin^ 2613 = 10"'^ and sin^ 2^12 = 0.8 we 
see that p ~ 72%, p ~ 28%. In effect the forward shock of such explosions acts like a window 
sweeping across either the neutrino or antineutrino spectra through which we observe an 
altered mixture of the emitted components. For neutrino or antineutrino energies below the 
critical energies this leads to an enhanced Ue or z/g flux as shown in Figs. fl22l) and fl23l) . 



C. A Heavy Water Detector 

The first detector type we consider is a heavy water, SNO like, detector and the signal we 
focus upon is the ratio of charged current (CC) - both from electron neutrinos and electron 
antineutrinos - to neutral current (NC) event rates. This ratio of event rates is the same 



signal considered by Schirato & Fuller [27 1 even though the "^He counters and then the NaCl 
inserted into SNO gave the detector some capability of distinguishing between the Ue and 
events. To some extent the uncertainty of the initial spectra can be removed by considering 
this ratio of event rates though by no means is it eliminated. For the results that follow we 



calculate the event rates using the Nakamura et al. [49(] neutrino-deuteron cross sections. 

The basic understanding of what happens to the v^. or v^. flux helps us understand the 
temporal evolution of the CC/NC ratio shown in Fig. (IMl) for the 2D simulation and 
for the ID simulations in Figs. ( l25i) . ( l26l) and ( 127|) . For all combinations of hierarchy, 
pinch parameter and explosion energy the ratio begins to vary at around 1 — 2.5 s into the 
simulation. The initial increase in the ratio is due to larger charged current rates (the neural 
current rate is flavor blind) caused by greater electron neutrino flavor content (for the NH) 
or electron antineutrino flavor (for the IH). The change is a clear signature of a transition 
from adiabaticity to non-adiabaticity that occurs when the forward shock reaches the H 
resonances for neutrinos with E ~ (9(10 MeV). As expected, the ratio varies more strongly 
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FIG. 24: The ratio of charged current event rates to neutral current event rates in a heavy water 
detector as a function of time for neutrino propagation through the 25° radial slice of the the 2D 
simulation. The top panel is for an inverted hierarchy , the lower panel for a normal hierarchy. In 
each the lines represent different spectral parameters: a = 1 (solid), a = 2 (dotted), a = 3 (long 
dashed), a = 4 (dash dot) and a = 5 (long dash-double dot). 
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FIG. 25: The ratio of charged current event rates to neutral current event rates in a heavy water 
detector as a function of time for neutrino propagation through the ID simulation with Q = 
1.66 X 10^^ erg. The top panel is for an inverted hierarchy, the lower panel for a normal hierarchy. 
In each the lines are: a = 1 (solid), a = 2 (dotted), a = 3 (long dashed), a = 4 (dash dot) and 
a = 5 (long dash-double dot). 
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FIG. 26: The ratio of charged current event rates to neutral current event rates in a heavy water 
detector as a function of time for neutrino propagation through the ID simulation with Q = 
3.36 X 10^^ erg. The top panel is for an inverted hierarchy, the lower panel for a normal hierarchy. 
In each the lines are: a = 1 (solid), a = 2 (dotted), a = 3 (long dashed), a = 4 (dash dot) and 
a = 5 (long dash-double dot). 
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FIG. 27: The ratio of charged current event rates to neutral current events rates in a heavy 
water detector as a function of time for neutrino propagation through the ID simulation with 
Q = 4.51 X 10^^ erg. The top panel is for an inverted hierarchy, the lower panel for a normal 
hierarchy. In each the lines are: a = 1 (solid), a = 2 (dotted), a = 3 (long dashed), a = 4 (dash 
dot) and a = 5 (long dash-double dot). 
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for an inverted hierarchy due to the greater change in the electron antineutrino spectrum 
when the profile becomes non-adiabatic than for the electron neutrinos as indicated by Figs. 
(122]) and (1231) . More detailed examination of the figures reveals that for the particular case 
of the weak explosion - Fig. fl2^ - there is evidence of the shock effects moving up through 
the neutrino spectra because at late times we see the ratio dips beneath it's t = 0.0 value. 

For the other simulations the CC/NC ratio exhibits a great deal of variability at the level 
of a few percent that continues for several seconds. We also notice that the variability of 
the Q = 3.36 x 10^^ erg and Q = 4.51 x 10^^ erg ID simulations ceased after ~ 5 — 6 s 
when the forward shock had swept through the spectrum and the reverse shock had turned 
around while the CC/NC ratio for the 2D simulation was sustained for the entire 10 s. While 
observing a period of fluctuations in the CC/NC ratio would indicate that the explosion and 
the density profile was dynamic it is difficult to be more quantitative about the presence and 
motion of the various features of the profile from the use of this ratio of event rates alone. 



D. Water Cerenkov Detectors 



The second example of a detector we examine is a water Cerenkov type since this is the 
dominant supernova neutrino detection technology. We focus upon the positron spectrum 
created by inverse jS reactions of antineutrinos upon protons. The temporal variations of the 
neutrino signal we have discussed in the paper would therefore require an inverted hierarchy 
in order to be visible in this signal. The positron spectrum seen in a detector, $e+(i?e+), is 
given by 

/7 
dE^.P.^J^ (23) 

where Np is the number of protons in the detector, F^^ is the electron antineutrino flux 
and da/dEf.+ is the differential cross section. For our analysis we used the cross sections of 
Strumia and Vissanti SOj. The positron spectrum allows us to see the effect of the density 
proflle upon the neutrino flux as a function of neutrino energy but, before we show our 
results, we flrstly consider the relationship between the neutrino energy and the positron 
energy and our ability to resolve the features of the neutrino spectrum. The positron energy 
and neutrino energy are closely correlated as shown in Fig. ( l28l) . Features in the neutrino 
spectrum at a particular energy will appear in the positron spectrum at more or less the 
same energy. But the positron spectrum is smoothed relative to the neutrino spectrum by 
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FIG. 28: The 'mean' positron energy (Eg) as a function of the neutrino energy (sohd hue) given 
by equation (16) in Strumia and Vissanti [50 ] . The dashed hue is {Eg) = Ey. 




FIG. 29: The Super-K energy resolution error, a, (sohd hue) and the range of positron energies, 
A, (dashed) as a function of the positron kinetic energy. 
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two effects: the distribution of positron energies for a given neutrino energy, and the energy 
resolution of the detector. For a given neutrino energy the positrons emerge with a range 
of energies which we may invert so that for a given positron energy, Ee+ , there is a range of 
neutrino energies. A, from which it could have arisen. This range^ is given by 



m:, — nip — ml + 2 nip ( \ 



2 2 

(24) 

2 \mp - Pe+ - Ee+ mp + Pe+ - Ee+ J 

where Pe+ is the momentum of the positron and m„, and m^. are the masses of the 
neutron, proton and electron respectively. A is shown in Fig. fl29|) as a function of the 
positron kinetic energy. And then for the energy resolution of the detector we used a filter 
so that the spectrum of positrons at energy Ee+ is given by 

^^EA ! dE[,^{E[,)W{E^AE'^,) 

and we take W{E^+\E'^+) to be a normalized Gaussian with a mean E'^+ and a variance 
cr^(E^+) given by 

/ 1 \ 0-3467 

0-319 -7^7-T7T - (26) 



E^+ — me ' \Ee+ — me[MeV] 
This formula is our fit to the energy resolution of the Super-K detector [53|]. The energy 
resolution of the detector, a, is also shown in figure fl29|) as a function of the positron kinetic 
energy. Both smearing effects, the range. A, of neutrino energies that can produce a positron 
with energy Ee+, and the detector resolution, a, increase with Ee+ and we notice that beyond 
Ee+ ~ 40 MeV it is the range A that determines that determines the scale of the features we 
can resolve. Below Ee+ ~ 40 MeV the detector resolution limits out ability to see features of 
the neutrino spectrum. Future detectors will improve upon the energy resolution used here 
and so - since this is also the energy range where we should focus our attention because it 
is where the electron antineutrino flux will vary the largest due to the changing adiabaticity 
- in all our figures we plot both the spectrum as it would be measured with the Super-K 
energy resolution and with a detector that has perfect resolution. Of course there is a third 
effect one needs to consider: the Poisson error due to the finite number of detected events. 
This does not smooth out features of the spectrum but rather make them difficult to detect. 
In what follows we shall ignore the Poisson statistics. 



The formula is valid for positron energies such that Pe+ + Ee+ < nip. 
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FIG. 30: The positron spectrum per unit tonne of detector as a function of energy at various 
snapshot times for the ID simulation with a total energy deposition oi Q = 1.66 x 10^^ erg. In the 
top panel the times are t = 0.0 s (solid) and t = 2.5 s (long dashed), in the middle panel t = 3.1 s 
(solid) and t = 4.0 s (long dashed), and in the bottom panel t = 5.0 s (solid) and t = 5.8 s (long 
dashed). In all panels the short dashed lines indicate the spectrum convolved with our adopted 
detector energy resolution. 
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FIG. 31: The positron spectrum in the case of the inverted hierarchy per unit tonne of detector as 
a function of time at various positron energies for the ID simulation with a total energy deposition 
of Q = 1.66 X 10^1 erg. The curves are ^ = 10 MeV (long dashed), E = 20 MeV (dash dot), and 
£^ = 30 MeV (soHd). The short dashed hues indicate the spectrum convolved with our adopted 
detector energy resolution. 

In Fig. fl30l) we plot the positron spectrum as a function of positron energy in the case of 
the inverted hierarchy at various snapshot times using the neutrino flux we calculated for the 
weakest ID simulation, Q = 1.66 x 10^^ erg. With this spectral information, we can begin 
to see the features of the explosion. The shock is clearly seen at ~ 10 MeV at t = 2.5 s in 
the perfectly resolved spectrum. By t = 3.1 s the sharp change in spectrum has moved up 
to ii^ ~ 25 MeV and by t = 4.0 s we begin to see the backside of the shock at ~ 10 MeV 
and again at t = 5.0 s at ~ 15 MeV and t = 5.8 s at E' ~ 20 MeV. The shock is a less 
dramatic feature in the spectra convolved with our adopted detector energy resolution but 
is visible particularly at t = 2.5 s. Given enough events even present detector technology 
would allow us to see the shock move up through the positron spectrum from which we 
could begin to glean information about the explosion. The explosion features also appear 
in the evolution of the positron spectrum at particular energies as a function of time. This 
is shown in figure fl3T|) where we plot the positron spectrum for E = 10 MeV, ii^ = 20 MeV 
and E = 30 MeV. The general trend for each energy is oc exp(— t/r) since this was the 
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time dependence we used for the emitted spectra but on top of this exponential decay the 
shock is obviously visible for all three chosen energies in both filtered and unfiltered curves. 
The width of the forward shock feature grows with E as we noted earlier. In Fig. (1321) we 
plot the positron spectrum for the more powerful explosion with Q = 3.36 x 10^^ erg at 
various snapshot times. The phase effects we saw in the crossing probability for this model 
appear also in the spectrum for perfect detector resolution. They are largest in the range of 
10 — 30 MeV where the difference between adiabatic and non-adiabatic electron antineutrino 
flux is largest and where the distribution of positron energies for a given neutrino energy 
is narrow. At higher energies the phase effects largely disappear because of the smoothing 
effect of the positron distribution. This smoothing is also seen in the plot of the spectrum 
at various energies as a function of time shown in Fig. fl5^ since the 20 MeV positrons has 
much greater variability than then 30 MeV positrons even in the filtered spectrum. We can 
even find the signature of a reverse shock that stalled and headed back to the core. In Fig. 
( !34|) we plot the positron spectrum for the more energetic Q = 4.51 x 10^^ erg ID explosion. 
The snapshot times are the same as those shown in Fig. ([8]) for the density and Fig. (|T9l) for 
the crossing probability for this model. The returning reverse shock appears at ~ 10 MeV 
at t = 4.5 s which matches what we expect from the crossing probability seen in the middle 
panel of Fig. (fT9|). 

Finally in Fig. fl35|l we plot the spectrum for the 2D simulation at various snapshot times 
selected so that they match as closely as possible the snapshots we have shown of the profile 
and/or the crossing probability. The forward shock and phase effects can be seen here too 
in the unfiltered spectrum and the inward moving reverse shock is responsible for the phase 
effects seen at ~ 15 MeV at t = 6.4 s and in Fig. (120|) . Other features upon the positron 
spectrum can also be related to various features shown in the crossing probability for this 
model. 

VI. SUMMARY AND CONCLUSIONS 

The next Galactic supernova has the potential to provide a great deal of information 
about both neutrino mixing and about the supernova. The temporal and spectral evolution 
of the signal is altered as the star explodes due to the changes in the density profile of the 
star. The forward shock should be easily detected even with present detector technology as 
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FIG. 32: The positron spectrum in the case of the inverted hierarchy per unit tonne of detector as 
a function of time at various positron energies for the ID simulation with a total energy deposition 
of Q = 3.36 X 10^^ erg. In the top panel the times are t = 1.5 s (solid) and t = 2.1 s (long dashed), 
in the middle panel t = 2.4 s (solid) and t = 3.1 s (long dashed), and in the bottom panel t = 3.6 s 
(solid) and t = 4.9 s (long dashed). In all panels the short dashed lines indicate the spectrum 
convolved with our adopted detector energy resolution. 
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FIG. 33: The positron spectrum in the case of the inverted hierarchy per unit tonne of detector as 
a function of time at various positron energies for the ID simulation with a total energy deposition 
of Q = 3.36 X 10^^ erg, and an inverted hierarchy. The curves are E = 20 MeV (dash dot), and 
E = 30 MeV (solid). The short dashed lines indicate the spectrum convolved with our adopted 
detector energy resolution. 

either a feature sweeping up through the positron spectrum in a water Cerenkov detector 
or as a sudden increase in event rates for a detector with no spectral information such as a 
heavy water detector. Once identified, and if the progenitor were known, the information 
can be used to measure the shock speed, the strength of the shock and potentially also 
the stalled shock time delay. The reverse shock can be spotted qualitatively as the cause of 
fluctuations in either event rates or upon the positron spectrum. There is also the possibility 
for the reverse shock to be seen if its motion is reversed and it heads back to the core. The 
neutrino signal may also contain hints that the supernova was aspherical as a qualitative 
change in the phase effects. 
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FIG. 34: The positron spectrum per unit tonne of water Cerenkov detector as a function of energy at 
various snapshot times for the ID simulation with a total energy deposition oi Q = 4.51 x 10^^ erg. 
In the top panel the time is t = 4.0 s, in the middle panel t = 4.5 s, and in the bottom panel 
t = 4.9 s. In all panels the short dashed lines indicate the spectrum convolved with our adopted 
detector energy resolution. 
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FIG. 35: The positron spectrum per unit tonne of water Cerenkov detector as a function of energy 
at various snapshot times for the 2D simulation. In the top panel the times are t = 1.6 s (solid) 
and t = 2.4 s (long dashed), in the middle panel t = 6.4 s (solid) and t = 7.4 s (long dashed), and 
in the bottom panel t = 8.0 s (solid) and t = 9.0 s (long dashed). In all panels the short dashed 
lines indicate the spectrum convolved with our adopted detector energy resolution. 
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